PROCEEDINGS 

OF SCIENCE 



Features of the Refined Gribov-Zwanziger theory: 
propagators, BRST soft symmetry breaking and 
glueball masses 



Silvio P. Sorella*+ 

Departamento de Fisica Tedrica, Instituto de Fisica, UERJ - Universidade do Estado do Rio de 
Janeiro, Rua Sao Francisco Xavier 524, 20550-013 Maracand, Rio de Janeiro, Brasil 
ii-OTfli'/; sorellaSuer j .br 

David Dudal 



(N 



Ghent University, Department of Physics and Astronomy, Krijgslaan 281 -S9, 9000 Gent, 
[t I Belgium 

£'-ma!7; david . dudal@ugent . be 

I — I Marcelo S. Guimaraes 

Departamento de Fisica Tedrica, Instituto de Fisica, UERJ - Universidade do Estado do Rio de 
JL Janeiro, Rua Sao Francisco Xavier 524, 20550-013 Maracana, Rio de Janeiro, Brasil 

^ Zi-zizfl/Z; msguimaraes@uer j .br 

I — I Nele Vandersicl<el 

^ Ghent University, Department of Physics and Astronomy, Krijgslaan 281-S9, 9000 Gent, 

^ Belgium 

\l E-mail: nele.vandersickel@uqent.be 

in 

The present work discusses an approach to access the physical spectrum of the Yang-Mills theory 
quantized in the Landau gauge. By using recent lattice data on the gluon propagator, it is possible 
to study the two-point functions of gauge invariant composite operators, from which masses of 
glueballs can be extracted. It turns out that the momentum dependence of the gluon propagator is 
very well reproduced by the corresponding tree-level gluon propagator obtained from the Refined 
Gribov-Zwanziger theory, which takes into account the presence of the Gribov horizon as well 
as the effect of condensates of mass dimension two. The resulting glueball masses are in good 



agreement with the available lattice data. 



The many faces ofQCD 
November 1-5, 2010 
Gent Belgium 



* Speaker. 

^Work supported by FAPERJ, Fundajao de Amparo a Pesquisa do Estado do Rio de Janeiro, under the program 
Cientista do Nosso Estado, E-26/101. 578/2010 



© Copyright owned by the author{s) under the terms of the Creative Commons Attribution-NonCommercial-ShureAHke Licence. 



http : / /pos . sissa . it / 



Features of the Gribov-Zwanziger theory and of its refined version 



Silvio P. Sorella 



1. Introduction 

The characterization of the physical spectrum of QCD at low energy scales from the first principles 
of the theory is an outstanding problem. QCD is a strongly interacting theory at low energy scales 
and thus the usual perturbative treatment does not apply. Due to color confinement, there is no 
free particle interpretation for the elementary fields appearing in the Lagrangian formulation of the 
theory. As a consequence, the properties of the physical spectrum are expected to be encoded in 
the correlation functions of suitable composite operators built with the elementary fields. 

The numerical lattice studies of the Yang-Mills sector of QCD, with the accurate determination 
of the ghost and gluon propagators, have provided valuable results. Gluons show unphysical prop- 
erties at low energies with a clear violation of positivity of the propagator. Although this behavior 
is expected in a confining theory, the employment of a positivity violating gluon propagator for a 
direct analytic evaluation of the correlation functions of the composite operators corresponding to 
the states of the spectrum is a highly nontrivial issue. 

In the following, we aim at showing that the decoupling type gluon propagator obtained from 
the Refined Gribov-Zwanziger theory yields a useful analytic set up for an estimate of the masses 
of the first glueball states with quantum numbers J^^ = 0++,0 ^,2++, while being in very good 
agreement with the recent lattice data. 

It is worth emphasizing that the framework we are going to outline here can, to some extent, be 
considered as model independent. In fact, the gluon propagator could be obtained by performing a 
good fit of the lattice data. The resulting expression could be employed to work out the correlation 
function of the composite operators corresponding to the aforementioned glueball states. Evidently, 
this is a model independent set up. It is remarkable thus that the gluon propagator obtained from 
the Refined Gribov-Zwanziger gives a very good fit of the lattice data in the non-perturbative low 
energy region, while providing a workable analytic structure. 

2. A short survey on the Gribov-Zwanziger theory 

The Yang-Mills theory is a gauge theory, thus, in order to be consistently quantized the gauge 
redundancy has to be removed. However, imposing a local, Lorentz covariant, condition which 
would select one and only one representative for each gauge orbit is not possible [1]. 

The problem was first noticed by Gribov [2] who observed that the Landau gauge condition d^A^^ = 
is not sufficient to account for the Gribov copies, i.e. gauge transformed configurations which 
still obey the Landau condition. Consider, for example, an infinitesimal gauge transformation 
— ^ +D"^(0^. The gauge transformed field would fulfill the Landau condition if d^D'^co'^ = 0. 
Zero modes of the Faddeev-operator 

^-1^ = -d,Df, Df = d,d'"-gf'"'A'^^ , (2.1) 

provide thus examples of Gribov copies. In order to take into account the existence of the gauge 
copies, Gribov [2] proposed to restrict the domain of integration in the functional integration to the 
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so called Gribov region Q.} : 

a=^Al\ d^Al = 0, > 0} . (2.2) 

The boundary of this region is known as the Gribov horizon and is the locus where the lowest 
eigenvalue of ^ vanishes. Later on, Zwanziger [3] was able to show that the restriction to the 
region Q. can be done within a local and renormalizable field theory framework. The resulting 
action is known as the Gribov-Zwanziger action 



Sgz = f d'xl^-F;i,F;i, + ib"d^A"^+c'^d^Dfc'^ 

+ Jd'x[fgf"^^A"^{(pf,^-^';)-d{N^-l)/) , (2.3) 

where d = 4 stands for the space-time dimensions. The dimensionful parameter 7 is known as the 
Gribov parameter. This expression can be written as 

Sgz = \J d^x f;^f;^ + sjd^x (c"d^A"^ - a)";dyDf(p'^''^ + Sy , (2.4) 
with Sy given by 

Sy = J d^x(^fgf"'A^;,{(p';-^';)-diN^-i)/) , (2.5) 

and where s denotes for the nilpotent BRST operator, whose action on the fields is specified by 



2 

sc" = lb" , 5^ = 0, 

s(pf = (Of!" , scof!" = . (2.6) 

The fields {(pf,'^f} and {(of,cof } are a set of bosonic, resp. fermionic fields introduced in 
order to express the action in local form. When the Gribov horizon is removed, which amounts 
to formally set 7 = 0, these fields give rise to a BRST quartet. They decouple from the theory, 
and the Gribov-Zwanziger action reduces to the ordinary Faddeev-Popov action. However, one 
has to observe that the parameter 7 is not free, being determined in a self consistent way by a gap 
equation, called the horizon condition [3]. It reads 

^ = 0, (2.7) 

'it is known that additional gauge copies still remain within the Gribov region. The elimination of all copies is, by 
definition, only attained by restricting the domain of integration to the Fundamental Modular Region. However, till now, 
a way to implement the restriction to the FMR within a local field theory framework is not available. 
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where Syac is the vacuum energy 

e-'~ = jme-'-, (2.8) 
and [(iO] stands for the functional integration over all fields appearing in Sgz- 

An important property of the action Sgz is its multiplicative renormalizability [3, 4, 5, 6] . This is 
a highly non-trivial feature that, among other things, enables us to consistently introduce into the 
theory composite operators and establish their renormalization properties. 

3. The issue of the BRST symmetry 

In the absence of the term Sj, the action (2.4) enjoys BRST invariance. In fact 

s j d'x Qiv^F;, + s (c-"5^A« - d)^' 5vDf ) = . (3.1) 

The Gribov-Zwanziger action is, however, not left invariant by the BRST transformations, eqs.(2.6), 
which are broken by the term Sy, namely 

sSgz = sSy = fJ d'x {-gf'"'iD"^^c')i(p'; - +gr''A^l,4') . (3.2) 

Though, the breaking of the BRST symmetry in the GZ framework is of a very special nature. 
Notice that the breaking term, being of dimension two in the fields, is a soft breaking. This fact en- 
sures the renormalizability of the theory through suitable Ward identities [3, 4, 5, 6]. Recently, there 
have been several developments in the understanding of this breaking. Notice that the breaking is 
quadratic in the fields. As such, it has to be treated as a composite field operator, a feature which 
requires the introduction of a suitable set of external sources in order to implement the Slavnov- 
Taylor identities. In [7], it has been shown that this breaking can be converted into a linear one. As 
a consequence, the resulting nilpotent Unearly broken BRST symmetry can be directly employed 
to obtain a system of Slavnov-Taylor identities. In [8, 9], it was pointed out that the softly broken 
BRST symmetry of the GZ action can be converted into a exact symmetry, however non-local. This 
non-local invariance has been localized in [10], though the resulting BRST symmetry is not nilpo- 
tent. Also, there are attempts to investigate the soft BRST breaking as a spontaneous symmetry 
breaking [11]. 

In summary, we can safely state that these results have already provided a certain understanding of 
the important issue of the BRST symmetry versus the Gribov horizon. The main result obtained 
so far remains that of the renormalizability of the theory. As a byproduct, let us mention the very 
important property that a local, gauge invariant operator ^ can be promoted to a renormalized 
operator ^r, so that the correlation functions 

(3.3) 

can be consistently evaluated order by order within the GZ framework [12]. 
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4. Dimension two condensates and the Refined Gribov-Zwanziger theory 

The action (2.4) is supposed to account for nonperturbative infrared features of Yang-Mills theory. 
The development of a non-zero value for the Gribov parameter 7 from the gap equation (2.7) is 
a manifestation of the nontrivial properties of the vacuum of the theory induced by restriction to 
the Gribov horizon. Moreover, there are additional sources on nonperturbative effects which can 
be encoded in a set of dimension 2 condensates, namely (A^Ap and ((p'^l'^'if — w^^ft)^''). It is 
possible to explicitly take into account the contribution of these condensates by introducing the 
con^esponding dimension two operators directly into the GZ action. The resulting formulation 
is, by now, known as the Refined Gribov-Zwanziger theory [6], and amounts to the following 
modification of the original Gribov-Zwanziger theory: 

Sgz^Srgz = Sgz + J d'xI^'^A^^A^^-M^ [r,'9f -^fcof)^ , (4.1) 

As well as the Gribov-Zwanziger action, the refined action Srgz also enjoys multiplicative renor- 
malizability [6]. 

5. A look at the gluon and ghost propagators 

Let us give a look at the gluon and ghost propagators obtained from both GZ and RGZ actions, 
respectively. Let us start with the GZ formulation. Here, the tree-level gluon propagator in <i = 4, 3 
and 2 dimensions is given by 

{A"^ik)A>^,{-k))Gz = 5"' (5,v - ^) , (5.1) 

where = Ig^Ny^. The gluon propagator turns out to be suppressed in the infrared region, attain- 
ing a vanishing value at zero momentum, /c = 0. On the other hand, the ghost propagator displays 
an enhanced behavior in the infrared 

(c-(%^-^))gz~5"^^, for ^2^0. (5.2) 

A gluon suppression and ghost enhancement are only observed for lattice data in 2 dimensions. 
They are not in agreement with the most recent lattice data for = 4,3 dimensions [13, 14, 15], 
which seem to point towards a suppressed gluon propagator which attains a non- vanishing value at 
k = 0, and to a non-enhanced ghost at ~ which keeps essentially the free behavior ~ ^ 

On the other hand, the gluon propagator obtained from the RGZ formulation in = 4,3 dimen- 
sions shows a different behavior. As in the GZ case, the propagator is suppressed in the infrared. 
However, differently from the GZ case, it attains now a nonvanishing value at A: = 0, as one observes 
from 

LAHt)4(-t)Uo. = - . (5-3) 
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where A"* = 2g^N'^ + m^M^. This expression is in very good agreement with the numerical data. 
In fact, ind = 4, an accurate fit with the SU(3) lattice data is possible for momentum scales up to 
A: ~ 1.5 GeV [16]. The fit is established for the following values of parameters 

M^ + m^ « 0.337 GeV^ « 2.15 GeV^ , A"^ ^ 0.26 GeV"^ . (5.4) 



Renorraalized Gluon Propagator - |4 = 3 GeV 
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Figure 1: Fit to the gluon propagator (left) and dressing function (right). Figures taken from [16]. 

Also, the ghost propagator obtained from RGZ in both d = 4,3 displays a behavior in harmony 
with the lattice data 

1 



{c^ikVi-k)) 



RGZ 



for k' ~ 0. 



(5.5) 



We also point out that in d = 2 the RGZ and the GZ theories coalesce, due to severe infrared 
divergences which do not allow to consistently introduce dimension two condensates. Due to this 
property, in d = 2 the RGZ theory gives the same results of the GZ theory. Therefore, the RGZ 
gives rise to gluon and ghost propagators which are in good agreement with the lattice data in 
(i = 4, 3, 2 dimensions. This is a remarkable achievement. 



6. /-particles and the glueball spectrum 

Both GZ and RGZ gluon propagators are associated to unphysical excitations. In fact, it is straight- 
forward to verify that the propagators (5. 1) and (5.3) each have two complex conjugated poles. Let 
us concentrate on the RGZ propagator. In this case, from (5.4), we can deduce the numerical values 
of the complex conjugated masses which, in GeV^ units, are 



:At^±\/202 = 0.1685 ±0.4812/. (6.1) 
The unphysical excitation associated with these complex masses have been called /-particles [17]. 



The /-particles diagonalize the quadratic part of the RGZ action, namely 

Srgz = J d'x ^-X; {-d^ + ml) X; + ^T]« (-^2 + ^2 ) r]^' +rest 



(6.2) 
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They provide a useful set up in order to extract the analytic properties of correlation functions of 
gauge invariant composite operators. For instance, the abelian part of the gauge invariant composite 
operator F^^F^y can be decomposed in the following way: 



If" f" ] 



= {d^Al-d,A"f = A^^T]" +rest. (6.3) 

abelian ft' 



where A^'y = dp^X" — dyX'^ and T]^^, = d^X^ — dyX'^ are the abelian field strengths of the /-particles 
fields. The first term, A^'vT]^^, of expression (6.3) can be shown to be associated to the part of 

the two-point coiTclation function of the operator {F"F") which displays a cut along the 

abelian 

negative real axis in the complex Euclidean plane, thus having a physical interpretation, see 
[17] for details. The remaining part, called rest in eq.(6.3), displays cuts along the imaginary axis 
and, as such, has no physical interpretation. Notice that the operator f^V^^'v has the quantum 
numbers of the glueball state J^^ = In order to study the glueball spectrum, in a lowest order 
approximation, we shall consider thus the abelian components of suitable composite operators 
^{x) displaying the quantum numbers J^'- and which, as in the case of {F"F") , are built 

I abelian 

with pairs of /-particles. Further, we look at the two-point functions of the operators i?{x). In 
order to have a physical meaning, these correlation functions should display a meaningful spectral 
representation, i.e. 

{&{k)i?{-k))= j^dtf^^. (6.4) 

with the spectral function p(?) ^ positive definite quantity. The spectral function has the general 
formp(f) =ZY,i5{t — nii) +Ad{t — to), whithZ andA positive quantities. The scales m, correspond 
to the poles of the correlation functions, thus identifying the masses of the physical particles of the 
spectrum. Moreover, is the threshold for a multi-particles state. We see thus that the spectral 
function carries important information about the physical spectrum of the theory. 
There are various methods to extract information on the spectrum from the knowledge of the spec- 
tral function, see, for example, [18]. In the present case, a modified reformulation of the moment 
problem [19], well adapted to the infrared region, was used to obtain estimates of the glueball 
masses of the scalar pseudoscalar ^ and tensor 2++ glueballs. The corresponding glueball 
operators in the /-particles representation are given by 

^o++(^) = K^{x)r]'l,y{x) + rest, (6.5) 



Pm«^v/3 - \PiivPa^ [Ka{^)^la{^)+^aa{^)^ai4) + rest, (6.6) 

^0-+ W = [^t^v{x)r]pa{x) + r]%{x)X"p^{x)j + rest. (6.7) 

Using the lattice input for the /-particle masses (6.1), the following estimates have been obtained: 
mo++« 1.96GeV, mo-+ 2.19 GeV , m2++ ?a 2.04 GeV. (6.8) 
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Note that the observed mass hierarchy mo++ < '«2++ < 'Mq y is reproduced^. By comparing these 

numbers with the lattice results 

rn^^l^ « 1 .73 GeV , m[fl+ « 2.59 GeV , mf_^^ w 2.40 GeV , (6.9) 

we can see a rather good agreement, all data being within a 20% range of approximation. 

7. Conclusions 

Our main conclusion is that the study of the glueball spectrum looks very promising within 
the RGZ framework. To some extent, our construction can be interpreted as showing a remarkable 
consistency of the most recent lattice data for the gluon propagator. These data have in fact been 
taken as input for the mass parameters entering the RGZ propagator which, in turn, has been em- 
ployed to achieve a rather good estimate for the masses of the three lightest glueballs. Also, the 
present results provide strong evidence that the interplay between the Gribov horizon and the di- 
mension two condensates, as described by the RGZ theory, accounts for a set on nontrivial aspects 
of the Yang-Mills theory at low energies. The RGZ gluon and ghost propagators turn out to be in 
agreement with the lattice data md = 4,3,2 dimensions. Also, these propagators have allowed us 
to obtain rather good estimates for the lightest glueballs starting from a first principle calculation 
of the correlation functions of suitable composite operators. 

Acknowledgements 

D.Dudal and N. Vandersickel are supported by the Research-Foundation Flanders (FWO Vlaan- 
deren). 

References 

[1] I. M. Singer, Commun. Math. Phys. 60 (1978) 7. 
[2] V. N. Gribov, Nucl. Phys. B 139 (1978) 1. 

[3] D. Zwanziger, Nucl. Phys. B 323 (1989) 513; ibid. 399 (1993) 477. 

[4] N. Maggiore, M. Schaden, Phys. Rev. D50, 6616-6625 (1994). [hep-th/9310111]. 

[5] J. A. Gracey, JHEP 0605, 052 (2006). [hep-ph/0605077]. 

[6] D. Dudal, J. A. Gracey, S. P Sorella et al, Phys. Rev. D78 (2008) 065047. [arXiv:0806.4348 
[hep-th]]. 

[7] M. A. L. Capri, A. J. Gomez, M. S. Guimaraes et al, Phys. Rev. D82, 105019 (2010). 
[arXiv: 1009.4135 [hep-th]]. 

[8] S. P Sorella, Phys. Rev. D80, 025013 (2009). [arXiv:0905.1010 [hep-th]]. 

[9] K. -I. Kondo, [arXiv:0905.1899 [hep-th]]. 

similar hierarchy has also been obtained in [20], where the ;-particles set up has been combined with the SVZ 
sum rules. 



8 



Features of the Gribov-Zwanziger theory and of its refined version 



Silvio P. Sorella 



[10] D. Dudal, N. Vandersickel, [arXiv: 1010.3927 [hep-th]]. 

[11] D. Zwanziger, Phys. Rev. D81, 125027 (2010). [arXiv: 1003.1080 [hep-ph]]. 

[12] D. Dudal, S. P Sorella, N. Vandersickel etal, JHEP 0908, 110 (2009). [arXiv:0906.4257 [hep-th]]. 

L13J A. Cucchieri, T. Mendes, PoS LAT2007, 297 (2007). [arXiv:0710.0412 [hep-lat]]. 

[14] A. Cucchieri, T. Mendes, Phys. Rev. Lett. 100, 241601 (2008). [arXiv:0712.3517 [hep-lat]]. 

[15] A. Cucchieri, T. Mendes, Phys. Rev D81, 016005 (2010). [arXiv:0904.4033 [hep-lat]]. 

[16] D. Dudal, O. Oliveira, N. Vandersickel, Phys. Rev D81, 074505 (2010). [arXiv: 1002.2374 [hep-lat]]. 

[17] L. Baulieu, D. Dudal, M. S. Guimaraes, M. Q. Huber, S. P. Sorella, N. Vandersickel and 
D. Zwanziger, Phys. Rev. D 82 (2010) 025021; arXiv: 1009.5846 [hep-th]]. 

[18] S. Narison, Camb. Monogr. Part. Phys. Nucl. Phys. Cosmol. 17 (2002) 1. 

[19] D. Dudal, M. S. Guimaraes, S. R Sorella, Phys. Rev Lett, to appear. |arXiv:1010.3638 [hep-th]]. 

[20] M. A. L. Capri, A. J. Gomez, M. S. Guimaraes etal, Eur. Phys. J. C 71 (2011) 1. [arXiv: 1009.3062 
[hep-th]]. 



9 



